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1 Introduction 



In the intermediate mass range of the Higgs-boson, its decay into photons is of great phenomenological 
interest. At hadron colhders the decay — > 77 provides precious informations for the discovery in the gluon- 
gluon production channel 1]. An upgrade option at the ILC will allow for a high precision measurement of 
the partial width into two photons ^ with a quantitative test for the existence of new charged particles. 

The QCD corrections to — > 77 have been computed in the past and analytic results at next-to-leading 
order are available in Ref. 3 and in Ref. [H (see also Ref. |5])- Electroweak two-loop corrections have been 
computed by suitable expansions of the two- loop Feynman diagrams 6 . Master integrals for the two- loop 
light fermion contributions have been analyzed in Ref. 7 and two-loop light fermion contributions to Higgs 
production and decays in Ref. 0. 

In our approach we have generated the full amplitude (up to two-loops and including QCD) in a completely 
independent way and we have used the techniques of Ref. [9] to produce a numerical evaluation of the partial 
width T{H — > 77). Since we are not bound to rely on expansion techniques, not even in the bosonic sector 
and in the top-bottom one, we can produce results with very high accuracy for any value of the Higgs-boson 
mass, taking into account the complete mass dependence of the -boson, Z -boson, Higgs-boson and top- 
quark. A consistent and gauge-invariant treatment of unstable particles made it also possible to produce 
very accurate results around the -threshold. 

2 Method of calculation and technical issues 

Our calculation builds upon the numerical approach of Ref. [9] where two-loop, two and three point 
functions have been investigated in the most general case. In this project we have developed a set of routines 
which go from standard Aq, ... , Dq functions up to diagrams needed for a two-loop 1 — > 2 process. This 
new ensemble of programs will succeed to the corresponding Library of TOPAZO jTU] . The whole collection of 
codes also uses the NAG- library [TT] . 

The generation as well as the manipulation of Feynman diagrams has been performed with the use of 
the FORM JJj code QraphShot [13]. Diagrams are generated, simplified and a FORTRAN interface is 
created. Furthermore, the code checks for the validity of the relevant Ward identities. Renormalization is 
performed according to the scheme developed in Ref. [13]. In this paper, we shall follow the same notations 
and conventions for two- loop diagrams as defined in Ref. |15j. In the following we will give a short outline 
of the techniques used for the calculation. 

Before evaluating the two- loop integrals arising after generating the Feynman diagrams, two main simpli- 
fications are done recursively. At first, reducible scalar products are removed and secondly, the symmetries 
of the diagrams are taken into account. The integrals are then assigned to scalar-, vector- and tensor-type 
integrals, according to the number of irreducible scalar products in the numerator and form-factors are intro- 
duced. The cancellation of scalar products is performed by expressing the scalar products in the numerator 
in terms of their associated propagators. This procedure can lead to removing lines in a diagram, so that 
each diagram produces a set of daughter-families with at least one line less. Apart from the reduction of 
scalar products, the consideration of the symmetries of a given diagram is important in order to reduce the 
number of integrals, which will be evaluated numerically at the end of the calculation. A simple example, 
showing the exploitation of the symmetries, is given in Fig. [1] for a scalar diagram. 

We now discuss briefly the extraction of coUinear logarithms from Feynman diagrams. It is worth noting 
that the amplitude for ^ 77 is collinear-free and one could adopt the approach where all light fermions 
are massless, then collinear behavior of single components is controlled in dimensional regularization and 
coUinear poles cancel in the total. We prefer another approach where collinear singularities are controlled 
by light fermion masses. Although the total amplitude is collinear-free, our procedure of reduction (g) 
symmetrization introduces a sum of several terms, of which some are divergent. Of course, we check that all 
logarithms of collinear origin cancel and, as a matter of fact, they cancel family by family of diagrams. 

To be precise, we need some universal representation for the coefficient of the collinear logarithms, which 
allows us to show their analytical cancellation, and a method to compute the remaining collinear-free parts. 
The first task is achieved by introducing integrals of one- loop functions and using their well-known properties 



P2 P'2 -P -P 




(a) Pi (b) Pi (c) Pi (d) 



Figure 1: Symmetries of the V^-family: The first diagram represents the 1^^-famiIy (a). Its integral remains 
unchanged by exchanging mi ^ m2 (b) as well as if one interchanges ms <-f m4 and p2 « ~P simultaneously (c). 
The last diagram (d) is a combination of the first (b) and the second (c) symmetry. One can also perform a total 
reflection of all external momenta, which is not shown in the figure and leaves the integral also unchanged. 



to make the cancellation explicit. Using the techniques of Ref. the coUinear finite contribution is first 
written in terms of smooth integrands and then evaluated numerically; an example is shown in Fig. [51 
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Figure 2: Example of a coUinear-divergent two-loop vertex diagram. Dashed lines represent particles with a small 
mass m and the wavy (external) line is massless. We introduced Lm ~ ln(m^/|P^|). 



3 Conceptual issues 

We will now apply our formalism to the computation of the amplitude for H{—P) + "f{pi) + 7(^2) — > 0, 
{P = pi + P2) which will be written as 
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16 7r2 



S^'' + 4^'^ P^ P", + 1^, Pl,P2) 



(1) 



The form factor is absent at O (g"^) and only arises at O (5^) but for a decay width with accuracy O (g^) 
(which includes one-loop (g) two-loop) its contribution is again zero. Bose symmetry and Ward identities 
(doubly-contracted, simply-contracted but with physical sources, simply-contracted with off-shell photons 
and unphysical sources) allow us to write the amplitude as 
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where the form factors are expanded up to two-loops, 
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A = A^^^ + 
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(21) 



(2) 



(3) 



2 



To proceed we need to include the relations between renormalized masses (small letters) and experimental, 
on-shell, ones (capital letters). Finite renormalization is then completed by introducing external wave- 
function factors {Zh^^'^ Z^^) and the renormalization of the coupling constants. All needed relations are 
collected in Eq.fg]). 
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where sl^l^, ^^hh and e|^^ are respectively the Higgs, W and top quark one- loop self-energies as defined in 



section 5.3 of the second paper of Ref. |T4]; furthermore, 
e1^|^.c(0) = £1^-1(0) + Jo, Sg^Q 
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The symbols Mt, M^, M^, Gp and a denote the mass of the top-quark, the -boson, the Z -boson as well 
as the Fermi-coupling constant and the fine structure constant. Collecting all the ingredients we get the 
corresponding S'-matrix completely written in terms of experimental data 
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The subscript "ex" indicates that all masses are the experimental ones and the mass-shell limit (s ^'l'^) 
is taken only after the inclusion of finite renormalization. QCD corrections will appear in Eq.(l4|) and Eq.® 
muhiplied by 7ra4M^)/(v^GpM^). 

In our calculation we prove the cancellation of the coUinear logarithms and then set the light fermion 
masses to zero; therefore, due to Yukawa couplings, an imaginary part in A^^^ arises only if M^ > 2M^. 
For two- loop terms imaginary parts are always present even for massless fermions. From Eq.(|6l) the total 
amplitude for iJ — > 77 can be written symbolically as -^p^ys ~ ■^1'^ <8i (1 + FR)+^2'^. Finite renormalization 
(FR) amounts to expressing renormalized parameters in the one-loop amplitude in terms of data and in the 
insertion of the Higgs wave-function factor Z^ a la LSZ; both requires the notion of on-shell mass. There 
are two sources of inconsistency in this approach: the Higgs-boson is an unstable particle and this fact has a 
consequence which shows up at two-loops. When we compute the doubly-contracted Ward identity for the 
full two-loop amplitude we obtain 



f almWiM„,M,,,...). 



(7) 



The analytical form of W is known and the non-zero result comes from the fact that the pure two-loop 
contribution to the Ward identity gives W while finite renormalization gives the real part Re W. Therefore, 
the Ward identity is violated above the VFVF -threshold. On top of this problem we find a second unphysical 
feature: let us analyze how the amplitude for iJ ^ 77 behaves around a normal-threshold, i.e. for M^ ~ 
2 M^,2 M^,2 Mt. In particular we are interested in the question of possible square- root or logarithmic 
singularities. Even if present they are unphysical, although integrable. Both problems can be solved by 
using complex masses as discussed in subsection [ 
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3.1 Square- root singularities 

It is very simple to prove that derivatives (represented by a dot in Eq.([H])) of one-loop, two-point functions 
with equal masses (to) develop a square root singularity: 

+ (8) 

The same argument can be repeated for all the one- and two-loop diagrams with any number of external legs 
where we can cut two and only two to- lines; normal-threshold will be a sub- (sub- . . . )leading singularity, but 
a 1/ /3-behavior shows up only if the reduced sub-graph responsible for the singularity can be reduced to a 
-Bq -function. Therefore the only two-loop vertex giving rise to a l//3-divergent behavior is the one depicted 
in Fig. [31 For this diagram it is possible to find a representation where the singular part is completely 
written in terms of one-loop diagrams, as shown in the figure. The remainder can be cast in a form suited 
for numerical integration. 



to 
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Figure 3: Contraction of a V^'' configuration leading to a /3 ^-behavior at the normal m-threshold. 

In the decay H 77, the 1 //3-singularity (/3^ — 1 — 4 M^/Af^ ) arises from the two-loop diagram of Fig. [31 
from Higgs-boson wave- function factor (derivative of a B-function) and from finite W -mass renormalization 
(derivative of a C-function). Our conclusion is that the unphysical l//3-behavior around some normal- 
threshold is induced by self-energy like insertion, a fact that is not surprising at all; those insertions, signaling 
the presence of an unstable particle, should not be there and complex poles should be used instead. 



3.2 Logarithmic singularities 

Let us consider the two-loop diagram of Fig. [H with = — s (s > 0). Writing the corresponding integral 



in parametric space we introduce the quadratic forms 

2 



4 TO 



and obtain 



=2 



^ dx dy 



Jo yx[x) 



Li2 1 



Li, 1 



yxjx) 
^(x,y) 



(9) 



(10) 



Since we are interested in the behavior around /3 — > 0, we split V'"^ into a singular and regular part and find 
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Figure 4: The irreducible, scalar, two-loop vertex diagram V"^ with logarithmic divergency. Solid lines represent a 
massive particle with mass m, whereas wavy lines correspond to massless particles. 



The singular part V^^^ will be written as [iGj 



V = — 

smg ^2 



dt- 



1-t 



lit), I{t) 



dx dy 



(1 - t) x{xy) + tyxix) 



dx dy 



a{x-xy + X 



where we have introduced the shorthands 
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with T = {\ — t)y + t and T ~ t/{l — t) > 0. I{t) can be split into two parts, 
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with Xi = —X, X2 = 1 — X and B(x,y) is the Euler beta-function. While the second term of Eq. fM]) is 
regular for (3 = 0, the singularity of the first term follows from the fact that A '--^ (1 — yY for (3^0; however, 
we have a singular behavior only if < X < 1 which requires y > ymm = max{0 , {t — l/2)/{t ~ 1)}. Since 
we are interested in the leading behavior for (3^0, we can extend the integration domain in the first term 
to [0, 1], without modifying the divergent behavior of the diagram. The singular part is then given by 
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< Res < 1/2. Here Fi denotes the first Appell- function. To obtain the expansion corresponding to (3 ^ 
we close the integration contour over the right-hand complex half-plane at infinity. The leading (double) 
pole is at s = 1/2. Therefore, we obtain 



J(t) = -- ln(-/32-iO)+C'(l), (3-^0. 
Inserting it into Eq.((n]) and using dtt^^^^ (1 - ty^^^ Int = - 27r, we get 

V^ns = ^ 1^ - ^ 0) + Oil), f3^0. 



(16) 



(17) 
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If the massive loop in Fig. 2] is made of top quarks the contribution of to the amphtude behaves hke 
/J^ and, therefore, the logarithmic singularity is /J^ -protected at threshold; the same is not true for a 
W-\oop. Our result, Eq.([T7]), is confirmed by the evaluation of of Ref. [T7] in terms of generalized log 
- sine functions. Starting from Eq.(6.34) of Ref. [IT] and using the results of Ref. [TH] we expand around 
9 = TT, where x ^ e^^ = (/? — 1) (/3 + 1), with < 9 < tt. This gives for the leading behavior of below 
threshold {tt'^ /i) \n{9 — tt), where ln(— /S^) = ln(9 — Tr)^ — In 2. The same behavior can also be extracted from 
the results of Ref. [4] . 



3.3 Complex masses 



Our pragmatical solution to the problems induced by unstable particles has been to remove the Re 
label in those terms that, coming from finite rcnormalization, give ReW^ in the Ward identity of Eq.([7|. 
Furthermore, we decompose Eq.((6]) according to: 
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and prove that, as expected. A},', A\' and AiJg satisfy (separately) the Ward identity. The latter fact 
allows us to - minimally - modify A^^^j^ by working in the complex-mass scheme of Ref. |19| . i.e. we include 
complex masses in the, gauge-invariant, leading part of the two-loop amplitude as well as in the one-loop 
part. 

The decomposition of Ea. (|18p deserves a further comment. There are three sources of 1//3 -terms: a) 
pure two-loop diagrams of the Vm -family, i.e. bubble insertions on the internal lines of the one-loop triangle; 
b) ly-mass renormalization, i.e. on-shell W-self-energy x the mass squared derivative of the one- loop W - 
triangle (the latter giving rise to 1//3); c) Higgs wave- function renormalization x lowest order (the former 
giving rise to 1//9). 

One can easily prove that only c) survives and a,b) that are separately singular add up to a finite 
contribution (/3 — > 0); their divergency is an artifact of expanding Dyson resummed propagators. 

The In [3 -term originates from pure two-loop diagrams (the Vk -family) and it is a remnant of the one-loop 
Coulomb singularity of one-loop sub-diagrams. 



4 Numerical results 

The partial width of the Higgs-boson decay into two photons can be written as 

The relative correction 5 induced by two-loop (NLO) effects is given by F = ro(l -f (5), where Fq is the lowest 
order result. It can be split into electroweak and QCF) contributions, 5 — f^Ew ^^qcd- 
For the numerical evaluation we use the following set of parameters: 

M„ = 80.398 GeV, = 91.1876 GeV, = 170.9 GeV, F,^ = 2.093 GeV, 

= 1.16637 X 10-5 GeV"^ a(0) = 1/137.0359911, as (M^) = 0.118. 

All light fermion masses are set to zero and we define the W -boson complex pole [20] by 

Sw = ^J-w ip-w - , fil, = M^^ -T'^, 7„, = Fv„- ( 1 - 2i|^) ■ (20) 

The one-loop H — > 77 amplitude, with a complex VF-mass, is shown in Fig. [5] around the M^VF -threshold 
including a comparison with the real W -mass amplitude. 
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Figure 5: Real and imaginary part of the one- loop // — > 77 amplitude with real and complex -boson mass. 
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Figure 6: Percentage electroweak corrections to H 
below the WVK -threshold. 



77 with real (dashed) and complex (solid) -boson mass, 



A comparison of the percentage electroweak corrections, with and without complex Vl^-masses, is shown 
in Fig. [6] for a Higgs mass range below the WW -threshold showing the unphysical growth of the real case 
and also some sizable difference in a region of about two GeV below the threshold. 

We have also analyzed the effect of (artificially) varying the imaginary part of the W -boson complex 
mass; results are given in Fig. [7l showing that our complex result reproduces the real one in the limit 0. 
Fig. [7] clearly demonstrates the large but artificial effects arising at normal-thresholds of unstable particles 
when their masses are kept real. 

Finally, in Fig. [8] we show both QCD and electroweak percentage corrections to the decay width T{H 
77), including the region around the PFW^ -threshold. A running as has been used for the computation of 
the QCD corrections. The remaining cusp of (5ew at the M^W^ -threshold, whose details are shown in the 

('2') 

blow-up of Fig. [U is due to our minimal scheme where the W -raaas is kept real in ^rcg, the regular part 
of the amplitude (see Eq. (fT8|) ). The relatively small error bars in a region so close to threshold serve as 
evidence for the efficiency of our numerical algorithms. 

Our result for Jew in the region 100 GeV < < 150 GeV is in substantial agreement with those of 
Ref. [6]. In conclusion, we observe a cancellation of the two corrections below the threshold whereas, above 
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Figure 7: The eflect of varying the -boson width in the real part of Aj^ = ^ ^ ~0~ containing the terms 
of the two- loop amplitude (left) and in the real part of V'^ of Fig. |4] containing the ln/3 terms (right) is shown. We 
also show the effect of using a real -boson mass but removing Re -labels in finite renormalization. 
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Figure 8: Electroweak (solid), QCD (dashed) and total (dashed-dotted) correction in percent for the partial width of 
the decay H —> 77. 



it, both (5qcd and S^w are positive leading to a sizable (up to 4.5%) total correction to the decay width. The 
perturbative expansion for the decay rate, supplemented with the complex-mass scheme, gives reliable and 
accurate predictions in a wide range of values for the Higgs-boson mass, typically —1% < Stot < 4% in the 
range 100 GeV < M„ < 170 GeV. 
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5 Conclusions 

In this paper we provide a stand-alone numerical calculation of the full two-loop corrections to the decay 
width T{H "f"f). Since no expansion is involved in the calculation we can produce results for all values of 
the Higgs-boson mass, as shown in Fig. [51 including the T4^M^ -threshold. The techniques introduced in this 
context are general enough to be used for all kinematical configurations of 1 ^ 2 processes at the two-loop 
level. 

To deal with normal-threshold singularities, specifically the M^VF -threshold, we have introduced complex 
W -masses in a gauge-invariant manner; our minimal scheme selects gauge-invariant components, typically 
LO (one-loop) amplitude and divergent parts of the NLO (two-loop) amplitude, and perform the replacement 
of Ea. (|^D)) . Details of our approach will be described in a forthcoming publication. 

The main result obtained in this paper can be summarized by saying that the NLO percentage corrections 
to the decay width T(H — > 77), (5qcd and (Sew, compensate below threshold leading to a small total correction; 
however, above the M^VF -threshold they are both positive, leading to a sizable overall effect of w 4%. 
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